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Abstract 

Large Mixing Angle (LMA) neutrino oscillation is the main solution for the long- 
standing Solar Neutrino Problem (SNP). Whether there is any subdominant effect 
accompanying the dominant LMA solution can not be ruled out at the moment, but 
will be settled by the forthcoming data from highly skilled real time experiments tar- 
geting essentially the low energy domain of solar neutrinos. Assuming a subdominant 
effect converting one of the active neutrinos into a sterile partner in the varying solar 
field with changing sunspot activity, we performed field-profile-independent predic- 
tions for ^Be neutrino signal variation, which might be tested in the KamLAND 's 
future solar neutrino detection program. We found that after a substantial reduction 
of background and running of KamLAND solar mode through the sunspot maximum 
period (around 2010 - 2012), when the solar field at the resonance may vary from few 
kG to 300 kG, the subdominant time variation effect might be clearly visible (more 
than 50") for ^Be neutrinos. 
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1 Introduction 



The KamLAND experiment [1] seems to confirm LMA as a solution to the decades long 
puzzle of missing solar neutrinos. Although LMA has been accepted as the dominant so- 
lution to the problem, there are several open questions which still remain to be answered 
[2]. Periodicity in the solar neutrino data, although controversial [3, 4], has been claimed 
recently by the Stanford Group [5, 6] upon examination of time binned data from all exper- 
iments except, so far, SNO. If it is confirmed, it will probably require Resonant Spin Flavor 
Precession (RSFP) within the sun affecting neutrinos in addition to the LMA effect. RSFP 
mechanism [7] converts active neutrinos of one fiavour into active or sterile (anti) neutrinos 
of another through the interaction of their magnetic moments with the transverse magnetic 
field on the way of their propagation. If the solar field is time varying this fiavour conversion 
is also time dependent. A significant time- modulation in the solar neutrino fiux would imply 
the existence of a sizable neutrino magnetic moment iii,, and hence a wealth of new physics. 

Most of the solar neutrino data analyzed till date [8, 9] is due to the high energy 
neutrinos (mostly E > 5 MeV), and in this way the low energy solar neutrino domain 
(essentially < 1 — 2MeV) is not much explored yet. In this sector the data we have so 
far are recorded by the inclusive radiochemical experiments Homestake, GALLEX/ GNO 
and SAGE. Examined in the time basis, some interesting signatures of time variability have 
been found in the data: The Stanford Group [5] claims a strong time variation effect for 
low energy neutrinos and, besides the unprovcn claim of time variation in Homestake CI 
experiment [10] ^, a time dependence in Ga data has been also noted [11]. The present 
statistics of experiments are however not enough to establish this effect. So, in order to 
confirm or rule out a time dependence in low energy solar neutrino sector we need real time 
experiments and enough statistics. 

Owing to the hints of flux variability signature suggested by the Gallium and Homestake 
Gl experiments, and as both the experiments have a sizable contribution of ^Be neutrino 
flux, it becomes worthwhile to study the ^Be neutrino line separately. Time variation for 
^Be neutrinos and other low energy neutrinos in the real time experiments BOREXINO and 
LENS has been examined in [12]. Our main motivation in this work is, as indicated by the 
title itself, to give predictions for the ^Be neutrino time variation in KamLAND experiment. 
KamLAND 's future plan is to detect and analyze solar active neutrinos in the low energy 
sector with ^Be neutrino measurements in the near future. 

The paper is organised as follows: in section 2 we discuss a possible time modulation of 
low energy solar neutrinos. The results and predictions for ^Be neutrino time modulation 
in KamLAND are exhibited in section 3, and flnally in section 4 the work is concluded. 

^The long time variation suggested several years ago for Homestake CI rate is still unclear [2]. 
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2 Low Energy Solar Neutrinos (LowE vs) 



Solar neutrinos arc overwhelmingly represented by the low energy neutrinos: 99% of total 
flux. The high energy tail (E > 5 MeV) of solar neutrinos contributes less than 0.01% 
of the total flux. LMA has been established essentially for high energy neutrinos because 
most of the solar neutrino data registered in the several detectors till date are due to these 
[8, 9]. On the other hand, LowE i/s, even though with dominating fraction, are almost 
concealed from the up to now studies of solar neutrinos. The quest for new physics in the 
LowE 1/ sector is at present the major objective for solar neutrino studies. 

2.1 Time Modulation in LowE us 

Apart from the unsettled claim of time variation in Homestake CI experiment [10, 2] Gal- 
lium experiments show a discrepancy of 2.4(T between the combined results of two different 
periods of time [11], which might be translated into a long time periodicity. Notice that Ga 
experiments have a dominant contribution of low energy neutrinos (pp Be ~ 80%). So, 
a strong time variation in Ga rates will definitely indicate time modulation in LowE i/ 
sector. 

In order to explain apparent signatures of time modulation in solar neutrino data [5, 6, 
10, 11] we developed a model [13], which predicts a clear time variational effect in the low 
energy sector of solar neutrinos and a hardly visible effect in the high energy tail of solar 
neutrino spectrum. In other words, model forebodes a maximum time variation of Ga rates, 
moderate of CI rates, and negligible effects on the SK and SNO rates. In the model, the 
magnitude of new mass squared difference (Am^Q) between one of the active states and the 
sterile one is chosen such that LowE us have spin fiavor precession resonating near the 
base of solar convective zone where the value of solar fields is peaked [14]. In this way there 
results a strong time modulation effect for these neutrinos. For qualitative illustrations see 
figure 1. To obtain this figure we used solar field profile 1 of [12] with an average field 
strength Bo = 150 kG. We plotted three probability curves showing strong suppression dips 
due to the subdominant RSFP mechanism in the low energy domain. The three dips (1,2,3) 
in the figure correspond to the three values of Am^Q, respectively as: 
1. -1.3 X 10-^ e\/2. 2. -1.0 X 10-^ eV^; 3. -0.6 x 10"^ eV^. 

Time dependence of LowE u fluxes is now an open possibility which might be settled 
by the transparent measurements with the forthcoming real time neutrino detectors like 
BOREXINO [15], LENS [16], and several others in the row [17]. KamLAND, a reactor 
antineutrino detector, is also programing to detect solar neutrinos, and would substantially 
contribute towards settling of this issue. 
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Figure 1: Time modulation effect in the low energy domain of solar neutrinos seen as dips 
(1,2,3) on survival probability curve due to the subdominant RSFP mechanism. 

3 ^ Be Neutrino Time Modulations in KamLAND 

KamLAND is aiming at the physics of LowE z/s, in particular "^Be neutrino hne, by 
lowering its threshold and eliminating the troubling background due to other sources. It 
will measure ^ Be signal via z/ — e elastic scattering (ES): 

u + e u + e 

where (') signifies the same particle but with energy changed. 

In order to perform predictions for ^Be neutrino time modulation in KamLAND, we 
estimated backgrounds and errors on the basis of informations collected from [18, 19]. For 
preventing the external gamma ray backgrounds from outside the liquid scintillator balloon 
we assumed its fiducial volume corresponding to a radius R < 4m leading to 7.2 x 10^^ target 
electrons. We assumed a substantial reduction of ^^Kr (/5-decays) background, although 
it is at present predominant in KamLAND over the Be neutrino signal region, and has 
to be reduced with a huge factor of order 10^. We also neglected the background due to 
^^°Po (a- decays), although it needs to be suppressed by a factor of ~ 10^, otherwise Be 



neutrino signal predictions become precarious essentially below 400 keV. We also assumed 
equilibrium after statistical subtraction of ^^^U, "^^"^Th backgrounds. By taking into account 
the smaller flux weights and negligible suppressions, we safely neglected the effect of time 
variations on the background due to pep, ^^O and ^^N neutrinos. This assumption is perfect 
for the energy window 400 keV to 800 keV. 

We used a Gaussian energy resolution function with cr{E) = Q.2%/\/E as published 
by KamLAND collaboration in [20], and adopted BP05(OS) [21] SSM flux value for ^Be 
neutrinos. Regarding the LMA parameters we considered KamLAND best fit for rate-and- 
shape analysis Ami^ = 7.9lg| x lO"^ eV^ and tan^0 = 0.46 [20]. We didn't use the combined 
data fit results of KamLAND + Solar Neutrino Data (Am|i = 7.9tg| x 10"^ eV^ and 
tan^O — 0.40lo;o7) because of the following reasons: 1) In the analysis of solar neutrino data 
only pure LMA has been considered and no subdominant RSFP effect is taken into account. 
2) The values of the parameters we used are within la of the combined KamLAND+Solar 
Analysis. Throughout the analysis we used fi^, = Ix 10^^^ fiB and allowed the average peak 
field value to vary in the range 10 — 300 kG. We chose the solar field profile 1 from [12], 
as preferred by the solar data fits, peaked near the bottom of convective zone as suggested 
by the observations [14]. 

3.1 ^Be Neutrino Spectrum and Model Predictions 

We studied the visible electron energy spectrum in the region from 280 keV to 800 keV, 
which is essentially dominated by the smeared ^Be neutrino signal, once the optimum re- 
duction of background is achieved. The visible electron kinetic spectrum is calculated as 

Res = Qo£'^ dT'fir, T)<PiE)[P,iE)^ + P^iE))^]- (1) 

Here Qo is the overall normalization factor, (p{E) is the SSM fiux of '^Be neutrinos and E is 
their energy. Pi{E) is the survival probability of fg, while P2{E) the production probabifity 
of Ufj, such that Pi{E) + P2{E) < 1 ^. The quantity f{T',T) is the energy resolution 
function of the detector in terms of the physical (T ) and the measured (T) electron kinetic 
energy, and and ^^^^ are the weak differential cross-sections for 1/^6 and f^e scatterings, 
respectively. 

Since ^Be neutrinos consist of a single energy line, by fixing Amfg all the neutrinos 
resonate exactly at the same position and experience the same field value. As a matter of 
fact, these neutrinos are blind to the shape of solar field profile, and in this way, because of 
this characteristic of '^Be neutrinos, the analysis is free from the uncertainties arising due to 

^Equality sign holds for pure LMA, when there is no time modulation. 
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the infinite possibilities of field profile shape^. In other words, the predictions for ''Be signal 
in KamLAND are quite robust in the sense that they are solar field-profile-independent. 

In order to work out the maximum possible time variation of ^ Be neutrinos we chose 
Am^Q = —1.3 X 10"^ so that ^ Be neutrinos resonate near the peak of solar field. At 
the sunspot maximum the peak field value could be as high as (300 — 500) kG [14]. We 
performed predictions for the electron spectrum (events/keV) in KamLAND for 3 years of its 
continuous operation, as these will preserve the long time variational effects due to sunspot 
activity, if any. In this way KamLAND could see a continuous increasing suppression up to 
the sunspot maximum (2010- 2012) and a continuous decrease in suppression subsequently. 
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Figure 2: Predicted visible electron energy spectrum events/ keV in KamLAND continuously 
running for 3 years. Lower curve refers to the spectrum due to the maximum possible time 
modulation and the shaded area around the curves refers to 99% CL. 

In figure 2 we have shown the expected spectrum in events per keV in the visible electron 
kinetic energy ranging from 280 keV to 800 keV . The maximum possible suppression is 

^The current knowledge about solar field profile is quite hazy; we have some observational support only 
for the peak field value Bo [14]. 



5 



shown by the lower curve with the average peak field Bo = 300 kG. The 99% CL errors are 
shown by the shaded region around the curves. The predictions of these curves and errors 
will be modified on the left of the vertical dotted line, once the ^^°Po background is taken 
into account. It is quite evident from the figure that in such a situation the time modulation 
effect during the solar cycle maximum (2010 - 2012) is clearly distinguishable (more than 
5a) in KamLAND data-taking for 3 years. For weaker fields {Bg ~ 140 kG) the two effects 
(pure LMA and LMA+RSFP) might be stiff distinguishable at 99% CL. 

After neglecting ^^°Po background and the time variation effects due to other solar 
neutrinos, which we did just for simplicity, the total number of events in the energy range 
280 keV to 800 keV expected in KamLAND 's three years of running is shown in table L 

We have chosen four test values for the average peak fields Bo- The first refers to pure LMA 
{Bo = kG), and the others to LMA + RSFP. In particular, the third value refers to a 
field of Bo = 140 kG (capable of separating time modulations from pure LMA oscillations 
at 99% CL), and the last value Bo = 300 kG refers to the maximum possible value of the 
peak field. 
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10^ 


LMA+RSFP 
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3.642 X 10^ 
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10^ 


3.642 X 10^ 


9.261 X 
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300 kG 


3.961 X 


10^ 


3.642 X 10"^ 


7.603 X 


10^ 



Table I - Total number of events in the energy range 280 keV to 800 keV expected in 3 
years of KamLAND data-taking. 

In order to avoid the contaminations due to ^^'^Po and the significant time variation 
effects on the background due to pep, ^^O and ^^N solar neutrinos, and to make our pre- 
dictions more tangible, it is legitimate to confine our analysis inside the window 400 keV 
- 800 keV. In this window KamLAND would only see the ^Pe neutrino signal once the 
backgrounds due to "^^^Bi and ^^Kr are significantly reduced. The total number of events 
expected in 3 years of KamLAND data-taking are shown below in table II. 
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400 keV 




800 keV 
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Bkgd. 


Total 


LMA 


kG 


4.285 X 10^ 


3.015 X 10^ 


7.301 X 10^ 


LMA+RSFP 


100 kG 
140 kG 
300 kG 


4.046 X 10^ 
3.836 X 10^ 
2.705 X 10^ 


3.015 X lO'^ 
3.015 X 10^ 
3.015 X 10^ 


7.061 X 10^ 
6.851 X 10^ 
5.720 X 10^ 



Table II - Total number of events in the energy range 400 keV to 800 keV expected in 
3 years of KamLAND data-taking. 




Figure 3: Individual ^Be and background visible electron energy spectrum events/ MeV/ sec 
in KamLAND. The dash-dotted spectrum refers to the maximum possible time modulation. 

The time modulation of individual ^Be neutrino spectrum in the energy range 280 keV 
to 800 keV is also shown in figure 3. The shaded area below the thin solid line as seen in 
the figure is the total background estimated in the KamLAND signal. The background is 
insensitive to the time modulation effect, which can be also concurred with the constant 
background event numbers given in the tables. However the ^Be neutrino signal is obviously 
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suppressed as the solar cycle evolves. In the figure the lowermost blue dot-dashed curve 
corresponds to the ^Be neutrino signal when = 300 kG (sunspot maximum), while the 
uppermost solid black corresponds to Bo = {0 — 10) kG (sunspot minimum, pure LMA). 
The blue dotted line refers to the suppression at 99% CL. 

It is evident from the figure that the lowermost curve is almost disappearing in the 
huge shade of background, implying that the KamLAND visible electron spectrum shape 
is being dominated by that of background, and in this way it would be quite hard to make 
out the shape ^Be neutrino spectrum. In other words, during sunspot minimum, which 
KamLAND might see around 2006 - 2007, the visible spectrum shape will be close to that 
of ^Be neutrino spectrum, however at the sunspot maximum the KamLAND integrated 
spectrum shape might be tampered by the background. In the later case if the background 
is significantly reduced, the dominance of '^Be neutrino spectrum shape could be retrieved. 
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Figure 4: Time variation in KamLAND with changing solar fields shown for three values of 
Amlo- 

In figure 4 we have depicted the total rate variation for three model values of Am^Q with the 
changing solar peak field value, with reference to pure LMA rate. We recall that the pure 
LMA solution corresponds to the field Bq — {0 — 10 kG) and, therefore, no modulations: 
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R/Rlma = 1, as seen in the figure. Modulations might push the pure LMA rate down 
by ~ 24% as shown by curve 1 exhibiting the maximum sensitivity to solar activity. In 
case of curve 2 the rate decreases slower, but faster than the curve 3 for increasing Bq. It 
is evident that curve 3 shows the least sensitivity to the solar field because ^Be neutrinos 
resonate far outside the base of solar convective zone and do not see an adequate field for 
their non-adiabatic conversion into sterile neutrinos. However, in this case pp neutrinos, 
having their resonance deeper inside sun around the peak field zone, may suffer maximum 
conversion into steriles, an effect which is beyond the current scope of KamLAND. 

4 Conclusions 

After an independent test confirmation of neutrino oscillation performed by KamLAND 
over a short basefine of ~ 200 kms, LMA has been accepted as a dominant solution to the 
SNP. Up to now most of the solar neutrino data are provided by the high precision real 
time experiments (SNO and SK), which have been monitoring the high energy tail of the 
solar neutrino spectrum. However, the dominant solar neutrino fiux (99%), which is due to 
low energy sector, is almost unexplored. On the other hand, the Stanford Group [5] claims 
even stronger time variation effects for LowE us. So, any subdominant effect in the low 
energy sector accompanying LMA can not be ruled out at the moment, because we don't 
have enough data and real time measurements which could establish or rule out any such 
effect. Nevertheless, such time variations, if they exist, could be certainly tested by the 
forthcoming high statistics real time LowE u experiments [17] . 

In this work we performed predictions for ^Be neutrino signal variation in KamLAND, 
assuming a possible correlation with sunspot activity. We assumed the ^Be neutrinos to 
resonate near the peak of solar field in order to examine the maximum possible time- variation 
in KamLAND. The results are quite robust for the reason that they are solar field-profile- 
independent. We have found that the time modulation effect might be clearly distinguishable 
from the pure LMA solution (more than 5a) in KamLAND data-taking for 3 years essentially 
corresponding to the sunspot maximum period around 2010 - 2012. It has also been seen 
that an average field of Bg — 140 kG can produce an effect capable of separating time 
modulations from pure LMA at 99% CL. At the time of low sunspot activity around 2006 
- 2008, the KamLAND electron spectrum shape is expected to be close to that of ^Be 
neutrino spectrum, and at the sunspot maximum the KamLAND integrated spectrum shape 
will be tampered by the background. Once the background is significantly reduced however, 
dominance of ^Be neutrino spectrum shape can be retrieved. Finally, we have also noted that 
at the time of sunspot maximum, the ^Be neutrino event rate might be further suppressed 
by ~ 24% relative to pure LMA. 
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